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1. Introduction 

Transport properties of degenerate dense matter have been studied for a long time. They 
are especially important for simulating various processes in neutron stars, white dwarfs 
and degenerate cores of giant stars. All these objects contain a degenerate electron gas, 
and the electrons can give considerable contribution to transport coefficients. Although 
the electron transport problem is very well elaborated and described in textbooks (e.g. 
PQ), some aspects have to be reconsidered. 

The first studies of electron transport properties (thermal conductivity) of 
degenerate stellar matter were performed by Marshak [2], Mestel [3] and Lee [I] in 
the 1940s. Further work in the next two decades is described by Lampe [5] who 
made considerable contribution into the calculation of the electron thermal conductivity. 
Later the subject has been studied by many authors, particularly, by Flowers and Itoh 
[6] and others O [HI [HI Ell ELT]. However, all these authors have considered collisions 
of relativistic charged particles in the static limit. The importance of dynamical 
interactions for relativistic particles was pointed out by Heiselberg and Pethick [T2] 
who analyzed transport properties of an ultrarelativistic quark plasma. The authors 
calculated transport coefficients of a degenerate ultrarelativistic plasma and showed that 
the results are qualitatively different from standard Fermi-liquid expressions because 
ultrarelativistic charged particles interact mainly through their currents. Such an 
interaction is produced via the exchange of transverse plasmons (instead of the standard 
Coulomb charge-charge interaction of non-relativistic particles via the exchange of 
longitudinal plasmons). An inclusion of transverse plasmons increases the effective 
collision rates (decreases the kinetic coefficients). 

Recently, we have reconsidered [13] the electron thermal conductivity of degenerate 
matter containing an electron gas of any degree of relativity. We have shown, that a 
correct treatment of electron-electron collisions, including the exchange of transverse 
plasmons, increases the contribution of these collisions into the electron thermal 
conductivity. 

In the present paper we perform similar revision of the electron shear viscosity. 
Such a viscosity is important for studying hydrodynamical processes in neutron stars 
and white dwarfs, particularly the damping of oscillations of these objects. We will focus 
on the contribution from the electron-electron collisions. Our analysis will be similar to 
that in Refs. [T2l fT3] . Thus, we omit technical details. 

The electron shear viscosity can be written as [1] 



where n e is the electron number density; p e = /i(37r 2 n e ) 1 / 3 and v e = p e /m* e are, 
respectively, the electron Fermi-momentum and Fermi- velocity, m* being the electron 
effective mass on the Fermi surface (it differs from bare electron mass due to relativistic 
effects). Furthermore, v e is the total electron effective collision frequency, a sum of 



n e v e p e 
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the electron-ion collision frequency is ei and the electron-electron collision frequency v ee . 
Detailed calculations of v ei have been recently performed in Ref. [13], while v ee is the 
main subject of our study. After calculating r] ee , we briefly analyze the efficiency of rj e in 
the cores of white dwarfs and envelopes of neutron stars (particularly, for the damping 
of pulsations of pre- white dwarf and white dwarf stars). 



2. Formalism 



We consider an almost ideal and uniform strongly degenerate electron gas. The electrons 
can have any degree of relativity and collide among themselves and with plasma ions. 
To calculate the electron-electron collision frequency we use the standard variational 
approach with the simplest trial function pQ. The variational expression for the effective 
electron-electron collision frequency, that determines the shear viscosity, is 

157r 2 h 3 f dpidp 2 dp' 1 dp' 2 



Uee 8v e pik B T J (2nh) 12 

x H/(12|l'2')/i/ 2 (l-/0(l-/2) 

X [pixVly + P2xV2y ~ Pl'xVl'y ~ P2'xV 2 'y} 2 , (3) 

where T is the temperature and k B is the Boltzmann constant. The integration is 
performed over all possible electron states involved into collisions P1P2 — > PiP 2 ; P is 
an electron momentum, p x being its component along the x axis, v y is the electron 
velocity component along the y axis; primes indicate particle states after a collision; 
/ is the Fermi-Dirac distribution function. Equation ([3]) includes the symmetry factor 
1/2 which excludes double counting of the same collision events of identical particles 
(electrons). Furthermore, W(12|1'2') is the differential transition probability, summed 
over spin states of colliding particles, 

x 5(p'i + p 2 - pi - p 2 ) i m /*i 2 ' ( 4 ) 

spins 

where e is the electron energy. The delta-functions reflect momentum and energy 
conservation; |M/j| 2 is the squared matrix element. 

The matrix element M/j for a collision of charged particles depends on the character 
of plasma screening of electromagnetic interaction between these particles. For ee 
collisions (of identical particles), one has M fi = MjJ } - Mf), where m}J } and ikfjf 
correspond to two collision channels, 1 — > 1'; 2 — > 2' and 1 — > 2'; 2 — >• V, respectively, 
and 

(i) _ 4?re 2 / jffiJyj _ Jfi'i ; Jt2'2 | /_v 

In this case, foq = p[ — pi and fko — e' x — S\ are momentum and energy transfers in 
a collision event, respectively; = (Jg°\ J e ' e ) = (2m*c)~ 1 (tte'7"tte) is the transition 
4-current, J te / e is the component of J e / e transverse to q; 7^ is a Dirac matrix; u e a 
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normalized electron bispinor (with u e u e = 2m e c 2 ), and u e is a Dirac conjugate. The 
expression for ML is obtained from Eq. (j3j) by the interchange of indices 1' <-» 2' (which 
results also in changing q — > q and uj — > —a;). 

For collisions of strongly degenerate particles, calculations are simplified by placing 
all interacting particles on their Fermi surfaces (whenever possible). Characteristic 
values of q and uj in electron-electron collisions are determined by plasma screening, 
that is described by polarization functions 11/ and H t for the longitudinal (charge- 
charge) and transverse (current- current) interactions, respectively. The nature of plasma 
screening is discussed in Refs. [121 13] • Let us summarize the main points. Usually, 
collisions of charged particles are studied in the so-called weak-screening approximation, 
in which momentum transfers are smaller than particle momenta. This approximation 
is justified by the long-range nature of electromagnetic interactions reflected in a specific 
g-dependence of the matrix element (jHJ) (a well pronounced peak at small q). As a result, 
only small values of q contribute to the integral ([3]), being determined by the character 
of plasma screening [12]. In the weak screening approximation, it is sufficient to consider 
the polarization functions in the classical limit (Hq <C p e and fru <C v e p e ), in which they 
are given by [15] 

n« = qlxi(x), IT = (q v e /c) 2 xt(x), (6) 
where x = uj/(qv e ), 

, . x 2 x(l - x 2 ) , fx + l\ ,„ N 
*<■*)= 2 +- L - r - Lin {j=l)' (7) 

h 2 qQ = 4e 2 pl/(irhv e ), and go is the Thomas- Fermi electron screening wavenumber. 
The longitudinal and transverse screenings are essentially different. The most striking 
difference occurs in a strongly degenerate plasma at temperatures T much below the 
electron plasma temperature T pe = hjj pe /k-B (that is determined by the electron plasma 
frequency u e = \jATte 2 n e /m* e ). In this case, it is sufficient to consider the low-frequency 
limit of uj —>■ and u/q v e , in which 

Xi = l> Xt = ircuj/(Aqv e ). (8) 

The longitudinal polarization function is real in this limit. It means, that longitudinal 
electromagnetic interaction (via the exchange of longitudinal plasmons) results in the 
standard Debye-like screening with a characteristic momentum transfer q\ ~ go [see 
Eq. (|5D]. For the transverse interaction (via the exchange of transverse plasmons), 
characteristic momentum transfer is different, q t ~ (tuq^/ve) 1 ^ 3 . Moreover, the transverse 
polarization function is pure imaginary. Accordingly, virtual transverse plasmons 
undergo collisionless absorption via the well-known Landau damping. 

Comparing screening lengths in the low-energy limit, we see that q t q\. Thus, 
transverse interactions occur on larger length-scales than longitudinal ones and are, 
therefore, more frequent. In the previous works on the electron shear viscosity, this 
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difference between longitudinal and transverse interactions have been neglected and one 
set II; = U t = q 2 . It is a good approximation in the non-relativistic case, because the 
ratio of transverse to longitudinal parts of the matrix element in Eq. (JSJ) contains a 
relativistic factor ( J te ' e / J^l) 2 ~ v 1l c2 - Hence, transverse interactions of non-relativistic 
particles are inefficient. In contrast, the collisions of relativistic electrons via the 
exchange of transverse plasmons are more important than the collisions via the exchange 
of longitudinal plasmons (owing to larger screening length). This effect was analyzed 
by Heiselberg and Pethick [12] for a gas of ultrarelativistic quarks. Here, we consider 
the electron-electron collisions for any degree of electron relativity, in analogy with the 
study of Ref. p3] for the thermal conductivity. 

Placing all colliding particles on their Fermi surfaces and performing possible 
analytical integrations with the aid of delta-functions, we finally obtain the following 
expressions for the electron-electron collision frequency and associated partial shear 
viscosity, 

12a 2 c 2 , 7ihn ePe v 3 e 

Here, a = e 2 /he is the fine structure constant, and 

f°° we w f 1 r def) 

I v (u,8) — / dw — / dx (1 — x 2 ) I — (1 — cos0) 



- I) 2 J, ' Jo vr 

1 u 2 (l — x 2 ) COS v 

(10) 



1 + (x6/w) 2 Xi( x ) 1 — u 2 x 2 + u 2 (x9/w) 2 Xt( x ) 
is a dimensionless function of two variables, 

u = v e /c, 9 = hv e q /(k B T) = \f?>T pe /T] (11) 

is the angle between components of momenta pu and p2t transverse to ^q, and w = 
huj / (ksT). Note, that in the small-momentum approximation, the interference term 
between and is small and both collision channels give equal contributions, 
resulting in \M fi \ 2 -> 2|m}J } | 2 . Equations © and $TD]) generalize Eqs. (47) and (48) 



from Ref. [12] to the case of arbitrary degree of relativity. 

Integration over in Eq. (11 01) is straightforward and gives 

I v (u, 6) = /,(«, 9) + I t (u, 9) + I ti (u, 9), (12) 

where // comes from the exchange of longitudinal plasmons, It from the exchange of 
transverse plasmons, and I t i is the interference term. 

Following Ref. [13] , consider four regimes of collisions between degenerate electrons 
as indicated in Table [TJ 

The regime I occurs in a non-relativistic (v e c) and rather warm (although 
degenerate) plasma (T > T pe ). The analysis of (fTUjl leads to the following asymptotic 
expressions valid in this regime: 

h =| Mn i + 1.919 
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Su 2 ( 1 

hi = ( In t: + 2.512 

15 V 

The leading contribution comes from // owing to a small relativistic factor u (as discussed 
above). The logarithmic terms and constant corrections in brackets are corresponding 
Coulomb logarithms (the corrections were calculated numerically). The leading-term 
result for r] ee in this regime is well-known. 

In the regime III, where the plasma is again warm [9 > 1), but the electron gas 
is ultrarelativistic (u ~ 1), Ii is the same as in the regime I (because Ii is independent 
of u). The asymptotic expressions for two other integrals have the same form as in the 
regime I but with different corrections, 

h=\ (ini + 2.742 

I tl = \ f In \ + 2.052 ) . (11) 



3 V 

Now all three terms give comparable contribution to I v . The asymptote for I v in the 
ultrarelativistic case coincides with that obtained by Heiselberg and Pethick [12]. It is 
different from the asymptote obtained by using an incorrect longitudinal screening in 
the transverse part of M/j. However, the difference occurs only in the corrections to the 
dominant logarithmic terms. 

The regimes II and IV are realized in a cold plasma < 1). In this case, we get 

r 3 



7T ' 



120' 

u 10/3 7T/4\ 1/3 



7T 



e-^T, Z=q(z) r(8/3)C(5/3) « 1.813 



7T 3 U 2 



- a - < 15 > 

where ((z) is the Riemann zeta function and T(z) is the gamma function. In the 
ultrarelativistic limit (the regime IV, u = 1), these expressions coincide with those 
obtained by Heiselberg and Pethick for an ultrarelativistic quark plasma [T2] . 



Table 1. Four regimes of shear viscosity i] ee of degenerate electrons. 



Electron Main T-dependence 



Regime velocity Temperature contribution of r/, 



I 


V e "C C 


T>T pe 


/i 


[Tln(T/T pe )] _1 


II 


V e C 


T<T pe 




1/T 2 


III 


V e ~ C 


T>T pe 


ii + + h 


[Tln^/T,,)]- 1 


IV 


f e ~ C 




h 


I/T 5 / 3 
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The main contribution to I v in the regime II comes from h due to a strong 
suppression of hi and J t in the non-relativistic case. In the regime IV, which operates 
usually everywhere in the neutron star crust (excluding only a thin layer - a few 
meters - from the surface) and in the neutron star core, the situation is different. The 
leading term is then I t , which corresponds to the exchange of transverse plasmons. It 
dominates because of smaller characteristic momentum transfers (larger electron mean- 
free paths) than those owing to the exchange of longitudinal plasmons. The temperature 
dependence of r] ee in the regime IV differs from the standard Fermi-liquid case (now 
i] ee oc T~ 5//3 instead of the standard dependence T -2 ). The difference is smaller, than 
for the thermal conductivity n ee (K ee is remarkably independent of T in the regime IV 
while a standard Fermi-liquid requires K ee oc T~ v ). Therefore, the Landau damping 
affects the shear viscosity weaker than the thermal conductivity. 

In addition to the above asymptotes, we have calculated J;, I t , and hi on a dense 
grid of u and 9 values covering transition regions between the regimes I-IV. We have 
further obtained the fit expressions which reproduce numerical results and asymptotic 
limits. The fit for h is 

0.287 \ . / 6.81' 



a379 + i + oW;o.ooi^ j'"( i+ ^j- < ie) 

with the maximum fit error of 1.4% at 9 = 1 (let us recall, that l\ is independent of u). 
The fit for I t reads 

1.813 (C 2 - 1.813/C)(1 + 53A9u - 764(0u) 2 ) 



C 



+ 



x In 



1 + 



1 + 7050m + 1630(#m) 2 
C 



(3 + 3.4m 2 )(#m) 



where G\ 



A(9uy/ 3 + (0u) 2 /3 
0.7801 + 0.1337m 4 , C 2 = 0.686 + 0.315m 4 , A = 0.556 - 0.08w 2 , 



(17) 
and 



C 



■ Aexp(Ci/C2)- The maximum error of 3.9% is at 9 
Finally, for I t i we obtain 



10" 



u 



0.57. 



hi 



5.168 (D 2 -5.168/D)(l+pi0 
D + ~ 



p 2 9 



2\1 



1 + Ps9 + p 4 9 2 + p 5 9 3 



x In 1 + 



D 
J 



0.354+1. 178m 2 , p 3 



where pi = 15.2 — 16.5m 2 , p 2 

p 5 = 0.00375 - 0.00373m 2 , D x = 1.339 + 0.029m 4 
D = exp(Dx/ D 2 ). The maximum error of 2.7% is at 9 



13.5-13.68m 2 , p 4 = 2.938-2.886m 2 , 
D 2 = 0.533 + 0.133m 4 , and 

10~ 2 , M = 1. 



3. Discussion 



Let us discuss the electron shear viscosity in a dense plasma. It is convenient to rewrite 
Eq. p as 



V, 

Vei 



-1 



= Vee + 

n e v e p e 
5z/« 



Vei ' Vee 



n e VePe 
5 V e .e. 



(19) 
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Figure 1. Temperature dependence of the electron shear viscosity at p = 10 6 g cm~ 3 . 
The line, marked 'ee', shows r\ ee \ the line 'ee-1' is the same, but only the contribution 
from the collisions via the exchange of longitudinal plasmons is retained. The dashed 
lines 'ei' and the solid lines 'tot' show r] e i and r] e , respectively, for helium and carbon 
plasmas. Vertical dotted line indicates the electron plasma temperature. 



where r] ee and r] e i are the partial viscosities governed by electron-electron and electron- 
ion collisions, respectively. We will discuss the effect of the Landau damping on rj ee and 
the contribution of r] ee in the total electron shear viscosity t] e . The partial viscosity T] e i 
will be calculated using the formalism of Ref. [13] . 

Figure [T] shows the temperature dependence of rj e in the helium and carbon plasma 
at p = 10 6 g cm -3 (where degenerate electrons become mildly relativistic). At higher 
p, the electrons are essentially relativistic and the effect of the Landau damping is most 
pronounced, see Section El Densities p > 10 6 g cm -3 are appropriate for degenerate 
cores of white dwarfs and red giants and for envelopes of neutron stars. We plot 
the partial shear viscosity rj ee (line 'ee') and the same viscosity, but retaining only 
collisions via the exchange of longitudinal plasmons (line 'ee-1'). These lines are the 
same for the helium and carbon plasmas. The relative contribution of collisions via 
the exchange of transverse plasmons in v ee increases when T falls below the electron 
plasma temperature T pe (log 10 T pe [K] rs 8.32, as indicated by the vertical dotted line in 
Figured]). The electron plasma temperature separates the high-temperature and low- 
temperature asymptotic regions (the regions III and IV for ultrarelativistic electrons, I 
and II for non- relativistic electrons, see Table [T]). 

The dashed lines in Figure [T] show i] ei and the solid lines marked 'tot' show the 
total electron viscosity r\ e . One can see that the electron-electron collisions are more 
efficient at higher temperatures in a plasma containing lighter nuclei (more exactly, the 
nuclei of lower charge, such as He, whose interaction with electrons is especially weak). 
At the given p = 10 6 g cm -3 , the plasma screening type is not very important, and the 
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Figure 2. Density dependence of the electron shear viscosity in a neutron star crust 
composed of the ground-state matter at T = 10 6 K. The line, marked 'ee', shows rj ee ; 
the line 'eeT is the same, but retaining the contribution from longitudinal plasmons 
alone. The dashed line 'ei' displays solid lines 'tot' give r/ e . The lines 'no US' and 'exp 
US' are calculated neglecting the Umklapp suppression and assuming an exponential 
Umklapp suppression in electron-ion collisions, respectively. The vertical dotted line 
shows neutron drip point. 

collisions via the exchange of transverse and longitudinal plasmons give comparable 
contribution. With increasing p, the electron-electron collisions become also more 
important at temperatures below T pe . However, at sufficiently high densities low- 
charge nuclei transform to those with higher charge (mainly owing to beta captures 
and nuclear fusion reactions). For these new nuclei with higher charge, the electron- 
electron collisions become less important. 

At high densities in a neutron star envelope (crust) the composition of matter 
cannot be arbitrary and is determined by an evolution scenario for a given star. One 
usually considers the models of ground-state (cold catalyzed) crust or accreted crust 
(e.g., Ref. [16]). Figure [2] shows the density dependence of r/ e , r/ ee and r] ei throughout a 
neutron star crust (for p > 10 9 g cm~ 3 ) at T = 10 6 K. We use a smooth composition 
model for the ground-state matter in the crust [IB] . Let us remind that the crust extends 
to p ~ 1.5 • 10 14 g cm -3 , while higher p correspond to a neutron star core composed of 
uniform neutron- rich nuclear matter. The vertical dotted line on the Figure [2] shows the 
neutron drip point (p n d = 4.3 • 10 11 g cm -3 ), which separates the outer crust (composed 
of electrons and nuclei) and the inner crust (where free neutrons appear in dense matter). 
The lines marked 'ee' and 'ee-1' again show r\ ee calculated including and excluding the 
exchange of transverse plasmons, respectively. Correct values of i] ee are more than one 
order of magnitude lower than the values 'ee-F. However, this suppression of i] ee by the 
Landau damping is insufficient for r] ee to dominate in rj e . The solid line marked 'no US' 
in Fig. [2] shows the total electron shear viscosity rj e in which r] ei is calculated under the 




10 r ' 

5 6 7 8 9 10 



log 10 T [K] 

Figure 3. Temperature dependence of the electron shear viscosity in a neutron star 
crust composed of ground-state matter at p = 10 13 g cm" 3 . Notations are the same 
as in Fig. [2] 

same assumptions as in [H]. Most importantly, the calculations neglect the freezing 
of Umklapp processes of electron-ion (electron-phonon) scattering at low temperatures 
(T <C T pe ) due to band-structure effects associated with motion of the electrons in 
crystalline lattice. It turns out that r] ei is several orders of magnitude lower than r] ee , 
and, therefore, r\ e rs r\ ei . 

However, the freezing of Umklapp electron-phonon scattering processes is a delicate 
task which has not been studied in detail in the literature. If the freezing operates, it 
enhances rj e i and makes rj ee more important. To illustrate this effect, in Fig.[2]we present 
the rj e and r] ei curves, marked by 'exp US' and calculated assuming an exponential 
freezing of the electron-phonon scattering rate after T falls below some temperature T u 
(that is approximately two orders of magnitude lower than ion plasma temperature). 
This illustrative model was used in Ref. [9] (and was described there in more detail). 
In this case the electron-electron collisions become significant at the bottom of a cold 
neutron star crust, at p > 10 13 g cm -3 . 

Figure [3] shows the temperature dependence of r] e , r] ei , and r] ee at p = 10 13 g cm -3 
(deep in the inner crust) under the same assumptions as in Fig. |2j If the freezing of 
Umklapp processes is neglected, then even for very low temperatures T ~ 10 5 -10 6 K 
(where r\ ee is actually two orders of magnitude lower than r] ee -i), i] ee gives negligible 
contribution to rj e (Figure [3j curves 'ee', 'ee-1' and 'tot', 'no US'). However, if one 
includes the freezing of Umklapp processes, then t] ee dominates at low T. 

The suppression of Umklapp processes at T < T u can actually be non-exponential, 
but power-law (A. I. Chugunov, private communication, 2007), which will reduce the 
importance of the electron-electron collisions. This suppression is complicated and 
should be a subject of separate study. In addition, we have neglected scattering 
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of electrons by charged impurities. It can also be important at low temperatures 
(depending on charges and abundance of impurity ions, e.g., [17]) and it can further 
reduce the importance of r\ ee (regardless the details of Umklapp freezing). 

The contribution of electron-electron collisions to r] e is not so high as their 
contribution to the electron thermal conductivity n e [13]. This is because K ee stronger 
depends on the screening momenta than rj ee . As a result, the dynamical Landau damping 
introduces into K ee an additional factor proportional to T. In the asymptotic region 
IV, K ee becomes then temperature-independent. In contrast, the Landau damping 
introduces into r] ee a much weaker factor oc T 1 / 3 . 

When calculating the electron-electron collision rate we have neglected the ion 
contribution into the polarization functions IT; and U t . It is a good approximation in 
the case of weak Coulomb coupling of ions T > Z 2 e 2 /(afcs), when the ions constitute a 
nearly ideal Boltzmann gas [a = (47mj/3) -1 / 3 being the ion-sphere radius, determined 
by the ion number density n^] . The calculation of the ion contribution to the plasma 
screening at lower temperatures T < Z 2 e 2 / '{aks) (where the ions form a Coulomb 
liquid or solid) is a complicated and unsolved problem. We have also neglected the 
effects of strong magnetic fields which can be available in neutron star envelopes and 
which can greatly modify the electron shear viscosity. These effects can be divided 
into two groups. First, they are the classical effects of electron magnetization owing 
to a fast electron rotation about magnetic field lines. Second, there are the effects 
of the Landau quantization of electron motion in a magnetic field (important usually 
for higher magnetic fields than the magnetization effects). The generalization of our 
solution to not too high magnetic fields, that do not affect the polarization functions, 
is straightforward. Stronger magnetic fields make the polarization tensor anisotropic, 
dependent of the magnetic field strength and direction, The effects of the ion polarization 
and strong magnetic fields on r] ee are beyond the scope of the present paper. 

4. Shear viscosity in the cores of pulsating pre-white dwarf and white 
dwarf stars 

Asteroseismology of white dwarfs is a rapidly developing field. More than 150 pulsating 
white dwarfs have already been observed (e.g., [18], [19] and references therein). The 
pulsation periods range from few minutes to few tens of minutes, and the relative 
pulsation amplitudes (pulsating fraction of star's luminosity) can be as high as a few 
percent. A comparison of observed and theoretical pulsation frequencies allows one to 
identify pulsation modes, to accurately determine white dwarf masses and radii, and to 
explore their rotation, magnetic fields, internal structure and evolution. 

All observed pulsation modes are interpreted as non-radial gravity modes (g-modes, 
produced owing to buoyancy forces) of multi-polarity £=1,2 with k = 1, . . .50 radial 
nodes. They are excited in the envelopes of young warm pre-white dwarfs and white 
dwarfs presumably by the instability that is mainly associated with partial ionization 
of the plasma. 




Figure 4. Kinematic shear viscosity rj/p in the core of the pre-white dwarf PG 1159- 
035 (left), the DBV white dwarf GD 358 (middle) and the DBA white dwarf R 548 
(right). The curves 'tot' show the total viscosity; the curves 'ei' and 'ee' are the 
contributions of electron-ion and electron-electron collisions, respectively. Lines 'ee- 
1' refer to electron-electron collisions retaining the exchange of longitudinal plasmons 
alone. The vertical scale on the right panel is different (shown on the right axis). 



Pulsating pre-white dwarfs (called PG 1159, or GW Vir stars) have the effective 
surface temperatures T e g from ~ 170 000 K to 75 000 K. They are hot, young stars 
(of age < 10 5 yr), still contracting slowly in the course of cooling because of a not too 
strong electron degeneracy in their cores. Their pulsations are thought to be driven by 
the partial ionization of C and O in the respective layers (which requires high T e g). In 
pulsating white dwarfs with helium atmospheres (DBV, or V777 Her stars), T e g- ranges 
from ~ 29 000 K to 22 000 K; these pulsations are excited in the He partial ionization 
zone. Pulsating white dwarfs with hydrogen atmospheres (DAV, or ZZ Ceti stars) have 
T e Q from ~ 12 500 K to 10 500 K. Their pulsations are most probably excited by 
convection in the outer layer with the partial ionization of hydrogen. Also, there exists 
a population of pulsating white dwarfs which belong to cataclysmic variables (accreting 
binaries). These white dwarfs have different composition of surface layers (contaminated 
by accretion from a companion star). As a result, they have different partial ionization 
zones, and their pulsations can be excited in a wide range of T e g (from 10 000 K to 
20 000 K and higher). 

Let us analyze the shear viscosity in the cores of pulsating white dwarfs. This 
viscosity participates in the damping of those pulsations which penetrate into the core. 
In a steadily pulsating star (whose pulsations are generated by a driving force) the 
viscosity can limit the pulsation amplitude. 

For illustration, in Figure H] we plot the density dependence of the kinematic shear 
viscosity rj/p in the cores of three pulsating stars (three panels) — one pre-white dwarf 
and two white dwarfs. A solid curve marked 'tot' on each panel presents the total 
electron viscosity rj/p, while dashed and dash-dotted curves 'ei' and 'ee' are the partial 
contributions of electron-ion and electron-electron collisions, respectively; a dash-dotted 
line 'ee-1' is the electron-electron contribution owing to the exchange of longitudinal 
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plasmons alone. 

The left panel of Figure H] corresponds to the conditions in the core of the 
extremely hot pulsating pre-white dwarf PG 1159-035. Its pulsation spectrum is 
rich; one has identified ^200 pulsation modes [20]. An analysis of astroseismological 
and spectral observations gives, among other things, the effective surface temperature 
T eff pa 140 000 K, and the star's mass pa 0.6 M [20]. According to the theory [21], the 
internal temperature of this star is mainly determined by the neutrino emission; it is 
nearly constant over the core, being close to the central temperature T c pa 10 8 K. Thus, 
we have plotted the shear viscosity for the isothermal degenerate core assuming pure 
carbon composition. The highest density p in the figure approximately corresponds to 
the central density of the star. At the lowest density, the electron degeneracy becomes 
mild (while our results are limited by a strong degeneracy). One can see that the 
electron-electron contribution to the electron shear viscosity is important only in the 
outer core, at p < 10 5 g cm -3 . The contribution of the electron-electron collisions due 
to the exchange of transverse phonons is small, but noticeable. It reduces the total 
electron viscosity r] e , maxiumum by a factor of 1.5, and the reduction increases with 
growing p. 

The middle panel of Figure H] is appropriate for the conditions in the core of the 
pulsating DBV white dwarf GD 358. It was the first DBV pulsating dwarf observed, 
with at least ten well-identified pulsation modes [22] (£ = 1, k < 20). Observations of 
GD 358 yield its mass ~ 0.6 M and the effective surface temperature T cS pa 24 000 K. 
These parameters and evolutionary models, reviewed in Ref. [23], imply the central 
temperature T c pa 4-10 7 K and the central density p c ~ 2-10 6 g cm -3 . We have evaluated 
the shear viscosity (Figure HJ) using an appropriate theoretical temperature profile in 
the white dwarf core from Ref. [23] . For illustration, we have selected the pure carbon 
composition of the degenerate core. One can see, that the electron-electron contribution 
r) ee to the total electron shear viscosity r] e is now less important, about 10 — 20%. The 
boundary of the degenerate core shifts now to lower densities (to p ~ 3 • 10 3 g cm -3 ); the 
contribution of the electron-electron collisions via the exchange of transverse plasmons 
is lower. 

Finally, the right panel of Figure [4] shows the shear viscosity for the conditions 
in the core of the pulsating DAV white dwarf R 548. The observational data imply 
T e g pa 12 000 K, and, again, the white dwarf mass ~ 0.6 M . The theory [23] predicts 
a nearly isothermal degenerate core with the temperature T pa T c pa 1.3 ■ 10 7 K. We 
have calculated the shear viscosity assuming the C-0 core with equal mass fractions of 
C and O. The shear viscosity r] e i, limited by electron-ion collisions in a C-0 mixture, 
has been determined using the linear mixing rule [16]. The contribution of collisions via 
the exchange of transverse plasmons is now noticeable only at p < 4 • 10 3 g cm -3 . 

Our three examples in Figure H] cover a representative range of pulsating pre-white 
dwarf and white dwarf models. In all the cases the kinematic electron shear viscosity 
r)/p is of the same order of magnitude, with characteristic values of rj/p ~ 10~ 3 — 10 -2 
cm 2 s -1 . If the damping of g- modes were solely determined by shear viscosity in the 
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stellar core, a typical damping time could have been estimated as r shear ~ X 2 p/rj, 
where A is a length scale of pulsation modes. Characteristic length scales are A ~ 
R/k ~ (0.1 — 0.01)i?, where R ~ 10 000 km is the white dwarf radius and k is the 
number of radial oscillation nodes. This gives r shear ~ 10 8 — 10 11 yr, indicating that the 
shear viscosity in the cores of pulsating white dwarfs is rather inefficient in damping 
observed pulsations. It would be equally inefficient to damp other possible large-scale 
hydrodynamic motions in white dwarf cores (for instance, differential rotation). 

The shear viscosity in white dwarf cores should have much stronger effect on small- 
scale motions. For instance, these could be small-scale (high-/c or high-£) pulsations, 
or ordinary pulsations in stratified cores, containing sharp boundaries between different 
layers (e.g., separating phases of different elements or phases of solidified and liquid 
matter). Viscous dissipation in boundary layers can be strong. 

Similar conclusions can be made on the efficiency of the electron shear viscosity 
in the envelope (crust) of neutron stars (where one typically has i] e /p ~ 0.01-100 cm 2 
s^ 1 for the range of the internal crust temperature from ~ 10 9 K to 10 7 K [Ej). The 
viscosity rj e there is not too high and cannot be a strong regulator of large-scale motions 
(pulsations). Nevertheless, characteristic length scales in neutron stars are naturally 
much shorter than in white dwarfs (because neutron stars are smaller, and their crust 
thickness is <1 km). This increases the efficiency of the shear viscosity (decreases 
viscous damping times) in neutron star envelopes (as compared to white dwarf cores). 
Moreover, a neutron star crust has a well defined heterogeneous (multi-layer) structure 
[T6] . where viscous boundary layers can occur and viscous dissipation can be strong. 
Such effects are currently almost unexplored. 

5. Conclusion 

We have calculated the partial electron shear viscosity r/ ee owing to collisions between 
degenerate electrons in a dense degenerate plasma taking into account the Landau 
damping of transverse plasmons. Our main conclusions are: 

(i) The Landau damping reduces i] ee for all temperature and density regions LTV 
(Tabled]). 

(ii) The strongest reduction occurs in the region IV of cold (T < T pe ) relativistic 
(p ^> 10 6 g cm -3 ) electron plasma. In this region, the Landau damping lowers 
the viscosity r] ee by several orders of magnitude and modifies its temperature 
dependence (which becomes p ee oc T _5//3 instead of the traditional dependence 
7] ee oc T~ 2 ). 

(iii) The viscosity r] ee gives a noticeable contribution to the total electron viscosity r\ e 
in a plasma of light ions at T > T pe . 

(iv) The viscosity r] ee can also give a significant contribution to r] e in a deep crust 
(p > 10 13 g cm -3 ) of a cold neutron star (T < 10 7 ) provided crystalline lattice 
of atomic nuclei is rather pure and the freezing of Umklapp processes at low T 
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is sufficiently strong. Both factors, the impurity of crystals and the character of 
Umklapp freezing, are currently uncertain and require special study. 

We have approximated our results by simple analytic expressions which are valid in 
the wide range of densities and temperatures appropriate to degenerate cores of white 
dwarfs and red giants and to the envelopes of neutron stars. The analytic approximations 
can be easily incorporated into computer codes to simulate hydrodynamical processes 
in these stars. We have briefly analyzed (Section H]) the efficiency of viscous damping in 
the cores of white dwarfs and envelopes of neutron stars. 

The contribution of collisions of degenerate electrons via the Landau damping 
during the exchange of transverse plasmons has been neglected in all previous 
considerations of the electron shear viscosity. In the context of transport properties 
of dense quark plasma it was studied by Heiselberg and Pethick [12]. We have recently 
reconsidered the effect of the Landau damping on the electron thermal conductivity 
in neutron star crust and core [TBI |2"1] . The same effect on the shear viscosity in the 
neutron star core (taking into account possible superfluidity of nucleons) will be studied 
in a similar way and published elsewhere. 
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